Mh in the MSSM-seesaw scenario with ILC precision 
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We review the computation of the one- loop radiative corrections from the neutrino/ 
sneutrino sector to the lightest Higgs boson mass, Mh, within the context of the so- 
i-C ' called MSSM-seesaw scenario. This model introduces right handed neutrinos and their 

supersymmetric partners, the sneutrinos, including Majorana mass terms. We find 
' negative and sizeable corrections to Mh, up to —5 GeV for a large Majorana scale, 

D , 10 13 — 10 15 GeV, and for the lightest neutrino mass in a range 0.1 — 1 eV. The corrections 

■ to Mh are substantially larger than the anticipated ILC precision for large regions of 

the MSSM-seesaw parameter space. 

1 Introduction 

The current experimental data on neutrino mass differences and neutrino mixing angles 
clearly indicate new physics beyond the so far successful Standard Model (SM) of Particle 
Physics. In particular, neutrino oscillations imply that at least two generations of neutrinos 
must be massive. Therefore, one needs to extend the SM to incorporate neutrino mass 
terms. 

We have explored [T] the simplest version of a Supersymmetric extension of the SM, the 
well known Minimal Supersymmetric Standard Model (MSSM), extended by right-handed 
Majorana neutrinos and where the seesaw mechanism of type I [5] is implemented to generate 
the small neutrino masses. For simplicity, as a first step, we focus here in the one generation 
case. 

On the other hand, it is well known that heavy Majorana neutrinos, with a Majorana 
mass scale to a/ ~ 10 13 — 10 15 GeV, induce large LFV rates [5], due to their potentially large 
Yukawa couplings to the Higgs sector. For the same reason, radiative corrections to Higgs 
boson masses due to such heavy Majorana neutrinos could also be relevant. Consequently, 
our study has been focused on the radiative corrections to the lightest MSSM CP-even 
h boson mass, Mh, due to the one- loop contributions from the neutrino/sneutrino sector 
within the MSSM-seesaw framework. 

In the following we briefly review the main relevant aspects of the calculation of the mass 
corrections and the numerical results. Further details can be found in pQ, where also an 
extensive list with references to previous works can be found. 
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2 The MSSM-seesaw model 



The MSSM-seesaw model with one neutrino/sneutrino generation is described in terms of 
the well known MSSM supcrpotcntial plus the new relevant terms given as: 



W 



Y V H\ UN -YiH\ UR + i N m M N, (1) 



where N = (v R , (vr) c ) is the additional superficld that contains the right-handed neutrino 
v R and its scalar partner v R . 

There are also new relevant terms in the soft SUSY breaking potential: 

Koft = m 2 L i>* L i> L + m R i> R P R + {Y v A v Hlv L v* R + m M B v v R v R + h.c.) . (2) 

After electro-weak (EW) symmetry breaking, the charged lepton and Dirac neutrino 
masses can be written as 

mi = Yi V\ , m D = Y y v 2 , (3) 

where are the vacuum expectation values (VEVs) of the neutral Higgs scalars, with 
^i(2) = v cos(sin)/3 and v ~ 174 GeV. 

The 2x2 neutrino mass matrix is given in terms of mo and mjvf by: 

M" = ( ° mD ) . (4) 

Diagonalization of M v leads to two mass eigenstates, which are Majorana fermions with the 
respective mass eigenvalues given by: 



m v , n = 7. ( m M T \l m 2 M + 4m i) ) ■ ( 5 ) 



In the seesaw limit, i.e. when £ = <C 1 . one finds, 



m v = -m D £_ + 0(m D £ 3 ) ~ , m N = m M + 0(m D £_) ~ m M ■ (6) 

rriM 

Regarding the sneutrino sector, the sneutrino mass matrices for the CP-even, M+, and 
the CP-odd, M_ , subsectors are given respectively by 

- 2 _ / m| + m 2 D + \M% cos 2/3 m D (A u - fxcot /3 ± m M ) \ , . 

± ~~ V m D{A u - [icotf3 ± m M ) m 2 R + m 2 D +m\ [ ±2B v mM ) 

The diagonalization of these two matrices, M±, leads to four sneutrino mass eigenstates. 
In the seesaw limit, where tum is much bigger than all the other scales the corresponding 
sneutrino masses are given by: 

m~ +fi _ = m| + -M\ cos2/3 ^2m D {A v -^ cot j3 - B v )£ , 
m% +! jf_ = m 2 M ± 2B v m M + m 2 R + 2m 2 D . (8) 
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Finally, in the interaction Lagrangian that is relevant for the present work, there are terms 
already present in the MSSM: the pure gauge interactions between the left-handed neutrinos 
and the Z boson, those between the 'left-handed' sneutrinos and the Higgs bosons, and those 
between the 'left-handed' sneutrinos and the Z bosons. In addition, in this MSSM-seesaw 
scenario, there are interactions driven by the neutrino Yukawa couplings (or equivalently 
tub since Y v = (gmo)/ (V^Mw sin (3)), as for instance ghv L v R = — % 2Mw sin p ' an< ^ new 
interactions due to the Majorana nature driven by mjf, which are not present in the case 
of Dirac fermions, as for instance ghv L v R = — 19 ™m™ sin °P " " Besides, the Higgs boson sector 
in the MSSM-seesaw model is as in the MSSM. 



3 Calculation 

In the Feynman diagrammatic (FD) approach the higher-order corrected CP-even Higgs 
boson masses in the MSSM, denoted here as Mh and Mh, are derived by finding the poles 
of the (h, _ff)-propagator matrix, which is equivalent to solving the following equation [4]: 

P 2 - m l + £ h h{p 2 ) 

where mh t H are the tree level masses. The one loop renormalized self-energies, E^(p ), in 
([9]) can be expressed in terms of the bare self-energies, £^(j> 2 ), the field renormalization 
constants 6Z$$ and the mass counter terms Sm 2 ^, where tfi stands for h, H. For example, the 
lightest Higgs boson renormalized self energy reads: 

t hh {p 2 ) = X hh (p 2 ) + 5Z hh (p 2 - ml) - 5m 2 h ■ (10) 

Regarding the renormalization prescription, we have used an on-shcll renormalization 
scheme for Mz,Myy and Ma mass counterterms and Th,Tn tadpole counterterms. On the 
other hand, we have used a modified DR scheme (mDR) for the renormalization of the wave 
function and tan /3 . The mDR scheme is very similar to the well known DR scheme but 
instead of subtracting the usual term proportional to A = - — 7s+log(47r) one subtracts the 
term proportional to A rn = A — \og(m 2 M / fjj£—) : hence, avoiding large logarithms of the large 
scale mjvf- As studied in other works [5], this scheme minimizes higher order corrections 
when two very different scales are involved in a calculation of radiative corrections. 

The full one-loop v/v corrections to the self-energies, E^", and entering ((9]) 

have been evaluated with FeynArts and FormCalc 0. The new Feynman rules for the 
vjv sector are inserted into a new model file. Since we are interested in exploring the 
relevance of the new radiative corrections to Mh from the neutrino/sneutrino sector, we will 
present here our results in terms of the mass difference with respect to the MSSM prediction. 
Consequently, we define, 

Am: M := M v h ' P - M h , (11) 

where M^ v denotes the pole for the light Higgs mass including the vjv corrections (i.e. in 
the MSSM-seesaw model), and Mh the corresponding pole in the MSSM, i.e without the 
vjv corrections. Thus, for a given set of input parameters we first calculate Mh in the 
MSSM with the help of FeynHiggs [7], such that all relevant known higher-order corrections 
are included. Then we add the new contributions from the neutrino/sneutrino sector and 
eventually compute Am" lDR . 



p 2 - m 2 H + Y, H h(p 2 ] 



t hH {p 2 ) 



= 



(9) 
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4 Results 



We have obtained the full analytical results for the renormalized Higgs boson self-energies 
and their expressions in the seesaw limit. In order to understand in simple terms the 
analytical behavior of our full numerical results we have expanded the renormalized self- 
energies in powers of the seesaw parameter £ = mo/wiM'- 



^-^U + ^U + ^U + ---- (12) 

The zeroth order of this expansion is precisely the pure gauge contribution and it does not 
depend on or mjvf. Therefore, it corresponds to the result in the MSSM. The rest of 
the terms of the expansion are the Yukawa contribution. The leading term of this Yukawa 
contribution is the 0{m 2 D ) term, because it is the only one not suppressed by the Majorana 
scale. In fact it goes as Y 2 M^ W , where Af| w denotes generically the electroweak scales 
involved, concretely, p 2 , M§ and M\. In particular, the 0(p 2 m 2 J ) terms of the renormalized 
sclf-cncrgy, which turn out to be the most relevant leading contributions, separated into the 
neutrino and sncutrino contributions, read: 



vmDR 
^hh 



-h + h--A V- 




g 2 p 2 m 2 D c 2 a g 2 p 2 m 2 D c 2 a 



(13) 

Notice that the above neutrino contributions come from the Yukawa interaction guu L v R , 
which is extremely suppressed in the Dirac case but can be large in the Majorana case. 
The sneutrino contributions come from the new couplings gho L c> R , which are not present in 
the Dirac case. It is also interesting to remark that these terms, being ~ p 2 , depend on 
the external momentum. Therefore, at large mj;, to keep just the Yukawa part is a good 
approximation, but to neglect the momentum dependence or to set the external momentum 
to zero are certainly not. In consequence, the effective potential method will not provide 
a realistic result for the radiative corrections to the Higgs mass. Similarly, obtaining the 
leading logarithmic terms in a RGE computation, would also miss these finite terms. 

The behaviour of the renormalized self-energy with all others parameters entering in the 
computation have been discussed in [T]. According to our detailed analysis in this paper, the 
most relevant parameters for our purposes are: mj; (or, equivalently, the heaviest physical 
Majorana neutrino mass tojv), Tfiv and the soft SUSY breaking parameters and B v . In 
the literature it is often assumed that rriM has a very large value, mj/ ~ 0(1O 14-15 ) GeV, in 
order to get \m v \ ~ 0.1 - 1 eV with large Yukawa couplings Y v ~ 0(1). This is an interesting 
possibility since it can lead to important phenomenological implications due to the large size 
of the radiative corrections driven by these large Y v . We have explored, however, not only 
these extreme values but the full range for mu'- ~ 10 2 — 10 15 GeV. 

Fig. [T] shows the predictions for Am" DR as a function of mu, for several input 
values. The Higgs mass corrections are positive and below 0.1 GeV if mu ^ 5 x 10 13 GeV 
and mjj < 10 12 GeV (left panel). For larger Majorana mass values, the corrections get 
negative and grow up to a few GeV; Am™ DR = —2.15 GeV for mu = 10 15 GeV. The 
results in the right plot show that for larger values of the soft mass, ms > 10 13 GeV, the 
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Figure 1: One- loop corrections to Mh from the neutrino/sneutrino sector as a function of 
m M for: < 10 13 GeV (left panel) and 10 13 GeV <m^< 10 14 GeV (right panel). 



Higgs mass corrections are negative and can be sizeable, a few tens of GeV, reaching their 
maximum values at ~ tum- 

The results of the Higgs mass corrections in terms of the two relevant physical Majorana 
neutrino masses, light \m v \ and heavy m^, are summarized in the left plot of Fig. [2] For 
values of < 3 x 10 13 GeV and \m v \ < 0.1 — 0.3 eV the corrections to M/j are positive 
and smaller than 0.1 GeV. In this region, the gauge contribution dominates. In fact, the 
wider black contour line with fixed Am™ DR = 0.09 coincides with the prediction for the 
case where just the gauge part in the self-energies have been included. This means that 
'the distance' of any other contour-line respect to this one represents the difference in the 
radiative corrections respect to the MSSM prediction. For larger values of uin and/or 
\m v \ the Yukawa part dominates, and the radiative corrections become negative and larger 
in absolute value, up to about —5 GeV in the right upper corner of this figure. These 
corrections grow in modulus proportionally to Tom and m v , due to the fact that the seesaw 
mechanism impose a relation between the three masses involved, m 2 D = \m v \rriN- 

Finally, we present in the right plot of Fig. [2] the contour-lines for fixed Am™ DR in the 
less conservative case where is close to ttim- These are displayed as a function of \m v \ 
and the ratio m^/rriM- m>M is fixed to 10 GeV. For the interval studied here, we see again 
that the radiative corrections can be negative and as large as tens of GeV in the upper right 
corner of the plot. For instance, Am™ DR = —30 GeV for \m v \ = 0.6 eV and m^/ruM = 0.7. 



To summarize: for some regions of the MSSM-seesaw parameter space, the corrections to 
Mh are of the order of several GeV. For all soft SUSY-breaking parameters at the TeV scale 
we find correction of up to —5 GeV to Mh- These corrections are substantially larger than 
the anticipated ILC precision of about 50 MeV (and also larger than the anticipated LHC 
precision of ~ 200 MeV). Consequently, they should be included in any phenomenological 
analysis of the Higgs sector in the MSSM-seesaw. 
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Figure 2: Contour- lines for the Higgs mass corrections from the v/v sector as a function of: 
|m„| and win (left panel) and the ratio m^/wiM and \m u \ (right panel). The other param- 
eters are: A v = B v = m L = m k = 10 3 GeV, tan f3 = 5, M A = 200 GeV, /i = 200 GeV. 
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